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Abstract. We present the results of our (the LU- 
NASKA project's) initial search for pulses of coher- 
ent radio Cherenkov radiation from UHE neutrino 
and cosmic ray interactions in the outer layers of the 
Moon. Three of the six 22 m antennas of the Australia 
Telescope Compact Array (ATCA) were implemented 
with our purpose-built pulse detection hardware 
and used to observe the Moon for two three-night 
periods during 2008. The observation dates and beam 
pointing position were chosen to maximise sensitivity 
to Centaurus A and the Galactic Centre (nominally 
Sgr A) region. We present the results of our (now 
complete) data analysis of order 10 7 candidate events 
- though we find no evidence for lunar Cherenkov 
pulses, our use of a wide (600 MHz) bandwidth has 
made our observations the most sensitive yet using 
the lunar technique. 
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I. Introduction 

The lunar Cherenkov technique fl] is a method to use 
radio-telescopes to detect ultra-high energy cosmic rays 
(CR) and neutrinos (v). By observing the short-duration 
(~few nanosecond) pulses of coherent Cherenkov radi- 
ation emitted from particle cascades via the Askaryan 
Effect [2 1 in the Moon's outer layers (nominally the 
regolith), the primary particles initiating the cascades 
may be identified. Our collaboration (LUNASKA) aims 
to develop the technique to be used with the next 
generation of giant radio-arrays. Here, we present the 
results of our two preliminary UHE particle searches 
using this technique with three antennas at the Australia 
Telescope Compact Array (ATCA) during February and 
May 2008. 

We leave to another contribution (Ref. [3 |) the details 
of our specialised equipment, and here concentrate on 
the results of our search for lunar Cherenkov pulses. 
We first describe the observations, in which we targeted 
Centaurus A and the Galactic Centre region, then detail 
the procedures used to discriminate between random 
noise events, RFI, and lunar-origin pulses. We calculate 
an effective aperture for our experiment to UHE v, and 
since our pulse search could identify no lunar pulses, this 
aperture is used to set an upper limit on an isotropic flux 



of UHE neutrinos. Finally, we calculate our exposure 
as a function of UHE v arrival direction, and compare 
this exposure with the exposure to Cen A and Sgr A 
from other experiments. Since we believe that current 
simulation techniques are inappropriate to modelling 
UHE CR hitting the Moon, we leave a calculation of 
our aperture/limit to UHE CR to a later work J4j. 

II. Observations 

We implemented three of the six 22 m ATCA antennas 
with our specialised pulse detection and dedispersion 
hardware, and observed over a continuous 600 MHz 
band in dual linear polarisations between 1.2 and 
1.8 GHz. This apparatus is described in greater detail in 
our other contribution (Ref. [3 1) - briefly, the equipment 
dedispersed, sampled (2.048 GHz, 8-bit), and scanned 
the data streams from both polarisations at each antenna 
for candidate pulses, and returned a data buffer of length 
256 samples upon detecting a signal greater than an 
adjustable threshold. A clock time accurate to 0.5 ns was 
also returned with each trigger. Each antenna operated 
independently due to data transmission constraints - 
since this limited our sensitivity, we set each antenna 
to trigger at approximately 100 Hz (5% dead-time) to 
partially compensate. 

A summary of our main observing runs and observa- 
tion time in February and May 2008 is given in Table 
|U The February run was tailored to 'target' a broad 
(> 20°) region of the sky near the Galactic Centre, 
harbouring the closest supermassive black hole to Earth 
and potential accelerator of UHE CR (or site of decaying 
super-heavy dark matter). For these runs we pointed 
the antennas towards the lunar centre, since this mode 
maximised coverage of the entire lunar limb, from which 
we expect to see the majority of pulses. In this mode, 
we achieved the greatest total effective aperture. 

Our May observing period targeted Centaurus A only, 
a nearby active galaxy which could potentially account 
for multiple UHE CR events as observed by the Pierre 
Auger observatory 0. Regardless of their source, this 
suggests the likelihood of an accompanying excess of 
UHE neutrinos, and we do not exclude the possibility 
of seeing the UHE CR themselves. We therefore pointed 
the antennas at that part of the lunar limb closest to Cen 
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Feb '08 1 


May '08 | 


Day 


26 th 


27 th 


28 th 


yjth 


18 th 


19 th 


t obs (mins) 


239 


319 


314 


324 


376 


440 


t efl p (mins) 


204 


277 


274 


224 


274 


316 



TABLE I 

Observation dates (nights thereof), total time t ohs , and 

EFFECTIVE OBSERVING TIME t eff AFTER SUBTRACTING TRIGGER 
DEAD-TIME, SPENT OBSERVING THE MOON IN DETECTION MODE 
OF OUR OBSERVING RUNS. 



A in order to maximise sensitivity to UHE particles from 
this region, as described in Ref. |4j. 

For each period, we observed nominally between the 
hours of 10 pm and 6 am local time. In 2008 we had to 
align the clocks manually by observing the bright quasar 
3C273 and correlating the emission between antennas, 
with a 0.25 ns timing uncertainty. Since T sys measure- 
ments were available over only part of our bandwidth, 
we used the Moon as our absolute sensitivity calibrator, 
with a 2% error. 

III. Data and Pulse search 

Our data set consisted of a few million triggers per 
antenna generated while pointing at the Moon. We there- 
fore wrote a simple algorithm to reject at least 99.99% 
of the random coincidences and RFI events, so the 
remaining data could then feasibly be inspected 'by-eye' . 
For this we used the timing criteria between antennas, 
since our three-antenna array allowed us to identify the 
arrival direction with extremely good accuracy. A pre- 
liminary search first identified all threefold coincidences 
- triggers on all three antennas within a time-envelope 
corresponding to the light travel-time between them - 
using the raw clock times at each antenna. Because 
of the large number of RFI events, this left of order 
100, 000 such coincidences. The relative clock times 
were then adjusted by pair-wise cross-correlations of 
the buffers, to allow for different antennas triggering on 
different parts of the signal. Events with a significant 
time-structure - e.g. Fig. Q] - produced a single strong 
correlation peak, with three pair-wise time adjustments 
that closed to zero. Those dominated by narrow-band 
RFI - e.g. Fig. |2] - tended to produce many peaks 
in the correlation. While this meant picking the right 
correlation peak was difficult, and usually the resulting 
adjustments in the offsets did not close to zero, this 
'problem' actually gave an easy method of eliminating 
such events. 

Using the adjusted time offsets and the known position 
of the Moon, a search over all remaining events produced 
60 coincidences where each pair-wise time offset was 
consistent with the lunar position, of which only 44 
resembled Fig. [2] and could be discarded. For our base- 
lines, the Moon was usually resolved by O ~ 30 'slices' 
perpendicular to the baseline, and of the remaining 16 
candidates, none produced three time offsets consistent 
with the same location on the Moon, and only one 




Fig. 1. Typical short-duration RFI event. 
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Fig. 2. Typical long-duration RFI event. 

was of sufficiently short duration to resemble a lunar 
Cherenkov pulse. Thus we concluded that no identifiable 
lunar Cherenkov pulses were observed. 

IV. Effective Experimental Aperture 

The instantaneous aperture, and consequent limits, 
from the experiment were calculated as per James and 
Protheroe 0). In all cases, we assume the existence 
of an effectively-infinite sub-regolith layer, and ignore 
secondary fi and r from the charged-current interactions 
of and v T . Since both the trigger threshold and differ- 
ence between assumed and actual STEC changed with 
time, so did the detection probability. We included the 
effects of finite sampling and approximate dedispersion 
by multiplying the simulated signal by a factor of order 
0.9, calculated as per Fig. 4 in Ref. [3|. Therefore, 
we calculated maximum and minimum apertures for 
each pointing configuration, and an approximate mean 
effective aperture, using the median values for the trigger 
thresholds. These are shown in Fig. [3] The detection 
threshold of a fewxlO 20 eV is comparable to that of 
the Parkes experiment, since the increased bandwidth 
compensated for the reduced collecting area. Also, it is 
lower for the May observations due to the lower T sys 
and greater sensitivity to the limb. However, at high 
energies, the ATCA observations (esp. for the centre- 
pointing configuration of February) see much more of 
the limb, and so are more sensitive than the Parkes 
experiment. 

Since we were targeting particular sources, the in- 
stantaneous 'beam-shape' (effective area as a function 
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Fig. 3. Isotropic apertures to UHE v, compared to the result for 
Parkes from [4], assuming the existence of a sub-regolith layer. 



of particle arrival direction) of the ATCA-Moon system 
was also of interest. Again using the methods of Ref. 
[4], we calculated this function for our observations, and 
plot it in Fig.|4] The limitation in targeting Cen A comes 
from the closest approach of the Moon being of order 
30°, where our sensitivity falls significantly, while for 
the Galactic Centre region, it is due to the Moon passing 
too close to Sgr A. Nonetheless, our exposure to these 
regions was significantly enhanced, as is shown in the 
next section. 

V. Limit on an Isotropic Flux 

Subtracting the trigger dead-time from the raw ob- 
servation times allowed us to calculate a 'model in- 
dependent' (2.3(A c st i )s )^ 1 ) limit on an isotropic flux 
of UHE v based on our null result. This is compared 
to existing limits in Fig. [5] The range for the ATCA 
observations reflects minimum and maximum sensitiv- 
ities; for the GLUE 0, Parkes Q, and Kalyazin O 
experiments (estimates from (4|), it reflects the inclu- 
sion/exclusion of a sub-regolith layer. While these limits 
exceed those of previous lunar Cherenkov experiments 
in the E v < 3 10 22 eV energy range (where NuMoon has 
stronger limits (9)), a 100-fold increase in observation 
time (approximately a full year's observation) would be 
needed at the current sensitivity to compete with the 
limits from ANITA iflOt 

While our isotropic limits are not competitive, the 
observation times and pointing positions were designed 
to target Cen A and the broad Sgr A region, as mentioned 
above. Integrating the instantaneous sensitivity of Fig. H] 
with time - given the known lunar position and orien- 
tation of the antenna beams - allows the total exposure 
(km 2 -days) as a function of celestial coordinates to be 
determined. This is plotted for neutrinos of 10 23 eV in 
Fig. [6j with approximate results for Parkes, GLUE, and 
RICE. The approximate declination range of ANITA is 
also shown, where the limit will dominate [10]. Table 
HIl compares the exposure from these observations to 
Cen A and Sgr A with that from RICE (Til, and our 
estimates for GLUE. Evidently, we have been successful 
in increasing our relative exposure to these potential 
sources of UHE particles. 




Fig. 4. Effective area (km 2 ) of the experiment to 10 22 eV neutrinos as 
a function of particle arrival direction, for (top) February, and (bottom) 
May. The telescope beams are centred on the Moon's centre (0, 0) and 
upper-right limb (0.17,0.17) respectively. 
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Fig. 5. Isotropic limits on an UHE v (see text). 



VI. Discussion and Conclusion 

We have successfully used an array of small, wide- 
bandwidth antennas to search for UHE particles using 
the lunar Cherenkov technique. That the limits on an 
isotropic flux are not competitive is expected, since 
these were only trial observations. However, the required 
observation time to set a competitive limit with the 
current system is probably impractical at this stage. 
We therefore concentrate our efforts on improving our 
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Fig. 6. Combined exposure (km 2 -days) from our LUNASKA ATCA observations, RICE | 111, Parkes |7*|, and GLUE (6) using our simulation 
results 0}, to 10 23 eV neutrinos. 
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E v 


GLUE 


RICE 


ATCA 


Total 


GLUE 


RICE 


ATCA 


Total 


10 21 eV 


0.5 


195 


2.9 


198 


0.015 


242 


6.9 


249 


10 22 eV 


14 


333 


54 


401 


2.1 


242 


111 


355 


10 23 eV 


175 


706 


409 


1290 


43 


512 


745 


1300 



TABLE II 

Experimental exposures (km 2 -days) of GLUE, RICE, and the LUNASKA ATCA observations to UHE neutrinos at 

DISCRETE ENERGIES FROM CEN A AND SGR A. 



sensitivity before conducting observations of serious 
length. 
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